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ANALYSIS FOR CONTROL APPLICATION CEF DYNAMIC 
CHARACTERISTICS OF TURBOJET ENGINE 
WITH TAIL-PIPE BURNING 
By Melvin S. Eeder and Richard Hood 

SUMMARY 

The general form of transfer functions for a turbojet engine 
with tail-pipe burning was developed and the relations among the 
variables in these functions were found from the transfer functions 
and from engine thermodynamics. By means of these relations, the 
dynamic characteristics of the engine can be found from steady- 
state data and one transient relation. 

The results of this analysis showed that if a step change in 
engine fuel flow causes the initial value of turbine -outlet tem- 
perature to be greater than the final value, a step change in 
exhaust-nozzle area or tail -pipe-burner fuel flow will cause the 
initial value of turbine-outlet temperature to be less than the 
final value, and conversely. Schedules that maintain constant 
engine speed and turbine-outlet temperature for a range of tail- 
pipe conditions and constant engine speed for a range of exhaust- 
nozzle areas were explicitly defined in terms of steady-state data. 

The results, when applied to the design of a noninteracting 
control system, gave the form of all the required control func- 
tions and showed that all but one of the control functions can be 
* determined from steady-state engine data. 


INTRODUCTION 

The application of tail-pipe burning to turbojet engines has 
become Increasingly important in improving the performance and the 
effectiveness of this type of aircraft power plant. The addition 
of the tail-pipe burner to the turbojet engine, however, increases 
the complexity of the control problem because of the additional 
degrees of freedom possessed by the engine and because both engine 
speed and temperature must be accurately controlled to obtain 
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niftTiTOnm engine performance with safe engine operation. Further- 
more, the action of one control variable during the transient state 
causes changes in the other variables. 

The first basic problem that must be solved before control 
synthesis can proceed is that of determining the dynamic character- 
istics of the engine. Accordingly, an an al ysis of the dynamic 
behavior of the turbojet engine- with tail-pipe burning was made at 
the BACA Lewis laboratory and is presented herein. 

In references 1 and 2, it is shown that the dominant dynamic 
characteristics of turbojet and turbine -propeller engines, respect- 
ively, may be expressed in terms of the slopes of engine-speed 
torque curves, which are derivable from engine performance data. 

It follows, as shown in references 1 and 2, that the engine proc- 
esses can be considered quasi -static, which implies that a thermo- 
dynamic process during transient conditions follows the path of 
equilibrium-state points. Thermodynamic relations that are valid 
for steady-state engine operation can therefore be extended to the 
transient state. Thermodynamic relations based on this result are 
used in reference 3 to determine the effect of the primary engine 
variables on the dynamic behavior of a turbojet engine with a 
centrifugal -flow compressor. 

In the present investigation of the dynamic behavior of a 
turbojet engine with tail -pipe burning, the general form of the 
engine transfer functions are developed and relations among the 
coefficients and the time constants are derived from the transfer . 
functions and from thermodynamic relations for the engine. The 
engine is considered a linear system in which incremental changes 
from steady-state operating conditions are considered. 

The relations among "the coefficients and the time constants 
are used to determine the indiclal response characteristics of the 
engine «nfl the analysis is then applied to scheduled and noninter- 
action controls . 


ANALYSIS 

General Form of Engine Dynamics 

The development and the data presented in reference 1 for a 
turbojet engine and in reference 2 for a turbine -propeller engine 
show that, at close-to-equilibrium operating conditions, unbalanced 
torque can be expressed as a function of tie engine speed and the 
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engine Independent variables . For a turbojet engine with tail-pipe 
burning operating at a constant ram pressure ratio and altitude, 
it therefore follows that unbalanced torque can be expressed as a 
function of the engine speed and the engine independent variables 
m in the following manner; (All symbols used in this report are de- 

o> fined in appendix A. ) 


Q = jr(*, F 0 , A, F t ) 


( 1 ) 


also 


Q <= I D(AE) (2) 

Equation (l), when expanded and linearized around steady-state 
operating points and combined with equation (2), leads to the 
transfer function for the response of engine speed to changes in 
the independent variables. This expression, which is developed in 
appendix B, is 


AN a l AF e | *2 AA **5 ^t 
N TD+1 F e TD+1 A TD+lF t 1 

As shown in appendix B, the response of turbine-outlet tem- 
perature to changes in the independent variables is given by 

,AT 2 T-jD+ 1 AF e TgD+l AA T 3 D + 1 AF t 

T 2 = TD+1 1 F 0 + TD+1 2 A + TD+1 3 F t 

The symbol A indicates incremental deviations from steady-state 
values. The variables F Q , A, and F^, which appear in the 
denominators, are the steady-state values. Thus, the term AN/E, 
for example, is the relative or percentage change in engine speed. 
The variables have been placed in this dimensionless form to make 
them consistent with the thermodynamic development to follow. 

From equilibrium conditions (D— > 0) and the principle of 
superposition, the coefficients a-j_ through b 3 are shown in 
appendix B to be proportional to the slopes of steady -state engine 
operating curves. These coefficients are defined in the following 
table ; 
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Coefficient 

Definition 

of 

coefficient 

Variables 

constant 

a l 

AN /AF e 
n/ F e 

A, F t 

®2 

AN /AA 
N/ A 

?e> ' 

a 3 

AN /^t 
N/F t 

*e> A 

*1 

AT 2 /"• 
W F e 

A, F t 

^2 

ATg /AA 

*r/ t 

F e , Ft 

*3 

AT 2 /AP t 
T 2 / »t 

F©> A 


The engine time constant t in equations (3) and (4) is a 
characteristic time in the transient solution of the homogeneous 
equations and can therefore he considered as the engine time con- 
stant -when the engine is displaced from, equilibrium and then re- 
leased (with all the Independent variables fixed). The significance 
of the time constants T l* T 2 > and T 3 may he illustrated hy di- 
viding equation (3) hy equation (4) after the response of If and Tg 
to the same forcing function is considered. For example, the response 
of N to Tg at constant A and F^ is 

ffl.-i-hiS. (s) 

If T-lD+ 1 h x Tg 

Thus, T]_ is the engine time constant at constant Tg, A, and 
Ffc. The engine time constants Tg and T 3 may he defined in a 
similar manner. The definitions of the engine time constants are 
presented in the following table: 
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Engine 

Variables 

time 

constant 

constant 


T 

F e , A, F t 

T 1 

^2* A> F-j. 

T 2 

V V p t 

T 3 

^ 2 * ^e* A 


Equations (3) and (4) present the general form of engine 
dynamics for a turbojet engine with tail-pipe "burning, "because, as 
shown in appendix B in the discussion of equation (B13), equa- 
tion (4) will be of the same form for any dependent variable. 

It also follows that the form of equations (3) and (4) is general 
for other engines in addition to turbojet engines with tail -pipe 
burning. For example, for the turbine -propeller engine, blade 
angle can be substituted for A in equations (3) and (4) and 
vanishes. The time constants and the coefficients oan then be 
defined in the manner previously described. 

A relation among the coefficients and the time constants in 
equations (3) and (4) is developed in appendix B. This relation 
is 

bi bp b* 

a ( V T) = A (t 2 -t) . 2 (Tj-t) (6) 

Because of the general form of equations (3) and (4), the relations 
in equation (6) hold for other engine typeB beside an engine with 
tail -pipe burning, a n d can be used by properly defining the time 
constants and the coefficients. 

Equations (3), (4), and (6) are developed without considering 
engine thermodynamics. In the analysis that follows, engine ther- 
modynamics will be used to develop relations among the time con- 
stants and the coefficients in equations ( 3 ) and (4), in addition 
to -Hie relation presented in equation (6). 


Engine Thermodynamics 

If the engine processes are considered quasi-static, thermo- 
dynamic relations that hold in the steady state can be extended to 
the transient state. Because unbalanced torque is a sma ll difference 
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'between two large members, only these thermodynamic express loos 
that are precise and Independent of engine efficiencies -will he 
used. These expressions are the heat -balance equation and the con- 
tinuity equation. 


Engine and tail-pipe equations . - For a turbojet engine, the 
heat-balance equation for the engine Is ■ 


He 



^ + H 2 

"a 


-Hi 


(7) 


where QU/tf a is the additional factor introduced by consideration 
of nonequilibrium conditions. 


The heat balance for a tail-pipe burner is 

Hfc ~ = Hg - Hg (8) 

"a 

The symbols Eq and H^ are introduced to make the equations 
that follow more compact. 


From the continuity of flow and the definition of Mach number, 
it follows that the gas flow through the exhaust nozzle is 


W, 


pgAM jy^s 

a “v^rvE 


( 9 ) 


Differentiation and linearization of equations . - If altitude 
and ram pressure ratio are assumed constant and if specific heat 
is assumed constant for differential changes in the variables, 
equations (7) to (9) can be differentiated to give the following 
expressions : 


Equation (7) becomes 




— dQ + 
F e 


H 2 dJ 2 
S e T 2 


and equation (8) becomes 

gt _ ga _ 

F t w a T 3 H*. T 2 


( 10 ) 


( 11 ) 
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* 


If the tail -pipe -"burner pressure ratio is assumed constant. 



( 12 ) 


Differentiation of the continuity equation yields 


dW a 

Wo 


cLA 


dP, 


A + P, 


1^ 
2 T« 




(13) 


A detailed development of equations (10) and (13) is presented in 
appendix B. 


The term dpg/pj in equation (13) can he omitted because, 
for constant altitude and subsonic flow in the exhaust nozzle, 
there is no change in static pressure, and for sonic flow the last 
term in the equation is zero. 

Equations (10) to (13) are in differential form. If it is 
assumed that the relations among the variables in equations (10) 
to (13) apply for incremental changes from steady -state conditions, 
the variables in these equations will be of the same linear form 
as those in equations (3) and (4). Thus, equation (10) can be 
considered as 


AF e AW a H t *2 AT 2 

F e "*a = F e + H e T 2 


(14) 


In this expression, AQ is the difference between the final and 
the initial unbalanced torque; therefore, because initial condi- 
tions are steady state, the initial unbalanced torque is zero and 
AQ. is equal to Q. In a similar manner, equations (11) to (13) 
nan be c ons idered in terms of incremental changes in the variables. 


Equations (11) to (13) can be combined to eliminate dP 3 and 
dT 3 . If these manipulations are performed and incremental changes 
from steady -state conditions are considered, the following equation 
results : 


5L\ 

AW a 

zeJ 

*a 


(1-hu) 


*2 



2H 3 E t 


(15) 
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where 



(The las t two terms in equation (15) sire bracketed hecause of ref- 
erence to them in the rest of the report.) 

In general, the term AWg/tfa in equations (14) and (15) is a 
function of engine speed and compressor-inlet a n d -outlet condi- 
tions. If an axial -flow compressor is assumed, however, the air 
flow may he approximated as a function of engine speed alone; 
therefore 


to 

a 





(16) 


u.Wq/ n« 

where W a ' = , which is proportional to the slope of the 

steady-state relation between engine speed and air flow. 

In the thermodynamic development, the variables have thus far 
been placed in the same form as those in the engine transfer func- 
tion (equations (3) and (4)). In order to obtain engine transfer 
functions expressed in operational form similar to equations (3) 
and (4), equations (14) to (16) are combined with one another and 
with equation (2) . 

Transfer functions . - Equations (14) to (16) cannot be used 
to completely describe the engine in terms of the independent vari- 
ables because APg.* which appears only in equation (15), cannot be 
eliminated from these expressions. A complete description of the 
engine would therefore require an independent expression for Pg. 

A physical consideration of the engine shows that P 2 can be 
expressed as a function of N, F e , and Tg. Because this func- 
tional relation does not contain A or F-^, and because the vari- 
ables AA arid AF-fc appear only in equation (15) and not in equa- 
tion (14) or (16), these variables will appear together only in 
the manner indicated in the bracketed term of equation (15). The 
engine time constants involving either of these variables are 
therefore equal, and relations among the coefficients in the trans- 
fer functions can be found in the following manner: 
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For convenience, equation (3) can fee rewritten as 

AN 1 . H TaA . f3 5£t] 

IT " TD+1 a l F e + TD+1 |_A + a F t J 

and from the discussion of equation (15), equation (17) 

M c 1 a ffe H f~AA _ JL. ^El\ 

IT C TD+1 1 F e TD+1 |_A " 2H 3 F t J 


(17) 


hecomes 

(Ifl) 


therefore 


f3 ^ _ ft_ 
®2 2E3 


Also, as shown in appendix B, the transfer function for the, 
response of IT to F e and Tg can he found by combining equa- 
tion (14) with equations (2) and (16) to eliminate AQ and AWa/ff a . 
This relation is 


AIT 1 ( **• AT 2^ 1 

N 0 Tj>D+l \F e "leT2/ W a * 


where, as shown in appendix B, 

t 2 bT 3 


ST 2 

»• V 


(19) 


( 20 ) 


Thus, from the thermodynamic development, the transfer function 
for the response of IT to changes in F e and Tj> is precisely 
defined in terms of steady -state variables. 


Equations (18) and. (19), when combined to eliminate 
give the following transfer function for the response of 
changes in the Independent variables: 




TD+1 


(l- a i V)He~ 
. H 2 


AF e 

»• 


ait/it, 
T 2 to 


■2 


V 1 
"a 


H e T 2 D+1 


AA 

A 


2H 3 Ffc 


H 2 TD+1 


(21) 
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find equation (4), rewritten in this form with t 2 equal to Tg, is 

( 22 ) 


ATg TjD+1 _ AF e 

(TgD + l\ 

AA ^t 

t 2 ~ TD+1 1 F e 

^tD+1 ) 

_ A + bg F t 


Equations (18) finfl (19), when combined to eliminate AF e /F e , 
give the following transfer function for the response of E to Tg, 
A, and Ft: 


AE 

E 


1 1 ) a l ^2 AT 2 + 

T ~*l V T g l-axW a *| H e T 2 + ^ 

V + 


AA 

A 


\ 


2H 3 F t 


(23) 


Equations (18), (19), (21), and (23) are transfer functions 
that describe the engine dynamics in terms of a minimum, of required 
data. 


Relations among coefficients and time constants . - Equa- 
tions (3) and (4) have been rewritten, as equations (17) and (22), 
to show relations among the coefficients and the time constants. 
Equations (21) and (22) are equal to one another and involve the 
same independent variables; corresponding terms are therefore 
equal. A similar correspondence exists between the terms of equa- 
tions (17) and (18). The following relations therefore exist 
between the time constants and the 'coefficients : 


_ T " a l W a* t 2 
1 l-a x W a « 




(1-ax W a ') H 


*e 
2 

b 3 „ , H e 

— = — = - Vi' — 

**3 ®2 ®2 




a 5 ^t 

ag b 2 2H 3 


J 


(24) 


and (repeated for convenience) the relations given in equations (6) 
and (20) 
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5l 

a i 


( t i- *0 


*2 

®2 


(t 2 -t) 



t 2 ssT 3 


in 2 

JeV 


( 6 ) 

( 20 ) 


Equations (6), (20), and (24) summarize the relations that 
have been developed among the time constants and the coefficients. 

Equation (6) : is not an Independent expression, hut can he 
developed from equations (20) and (24) and is consistent -with them. 
Equation (6) Is presented because, as indicated in the discussion 
of the equation, the relation is of a general nature. 


RESULTS AND DISCUSSION 


Engine Dynamic Characteristics 

The transfer functions of the engine are of a form for which' 
specific dynamic characteristics can he obtained. One such dynamic 
characteristic that is useful in controls analysis is the response 
of a system to a step input (indlcial response). Another import- 
ant characteristic is frequency response. 


The transfer functions of the engine are of two general forms, 
which are 


AY = oD+1 AX 
Y * pD+1 ° X 


(25) 


and, if a is zero, 


AY = 1 0 AX 

Y “ PD+1 X 


(26) 


Inasmuch as the properties of these transfer functions are well 
known, only the response of the functions to step changes in the 
Independent variables will he discussed in detail. 


Indlcial response . - For step changes in the independent var- 
iables, transfer functions of the form of equations (25) and (26) 
lead to exponential curves similar to those presented in figure 1. 
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Final equilibrium conditions for a step Input In the Independ- 
ent variables con be found from equations (25)' and (26) by allow- 
ing D to approach zero and are c as shown In figure 1. 

Initial conditions for a step Input can be found from these equa- 
tions by allowing D to approach Thus, for transfer func- 
tions of the form of equation (26), the indie ial response will be 
as shown in figure 1(a), where the, initial response to a step 
input is zero. If a/p >1 in the transfer functions of the form 
of equation (25), the initial value of the dependent variable will 
be greater than the final value and the lndicial response will be 
as shown in figure 1(b) . Similarly, if a/p < 1, the initial 
value of the dependent variable will be less than the final value 
and the lndicial response will be as shown in figure 1(c). 

The coefficient c in the transfer functions is a scale fac- 
tor that could be incorporated in the ordinates of the curves of 
figure 1. In the general case, the coefficient o can be posi- 
tive or negative and the ordinate can be interpreted consistent 
with the sign on this coefficient. 

The transfer functions of the engine, as derived in the pre- 
vious section, can now be interpreted in terms of the curves of 
figure 1. From equations (3) and (4), the Initial and final val- 
ues of the dependent variables for step changes in the independent 
variables can be found in the manner discussed for equations (25) 
and (26) . The results of this procedure are presented in the fol- 
lowing table; 


Step change 
in 

Variables 

constant 

Response of 
engine speed 

Response of 
turbine-outlet 
temperature 



Initial 

value 

Final 

value 

Initial 

value 

Final 

value 

AF e 

*e 

A, Ft 

0 

AP e 

ai i r 


AF e 

AA 

A 

Ft 

0 

AA 

•* T 

T 2 . AA 

T 

AA 

b 2T 

AF t 

Ft 

3*6, A 

0 

AF t 
®3 p t 

t 3 . 

T *3 , t 

AP t 

b3 pT 


This table shows that the lndicial response of engine speed 
to changes in the independent variables will be of the form given 
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in figure 1(a) . The indicial response of turbine -outlet tempera- 
ture to changes in the independent variables will he of the form 
shown in figure 1 (h) or 1 (c), depending upon whether the ratio of 
the time constants is greater or less than 1 . 


In figure 1, no consideration was given to the possibility of 
negative values of a. For a to.be negative, Tq, Tg, or T 3 
must he negative. From equation (20), however, Tg or T 3 will 
not he negative unless W a * is negative; ¥ a ’ is proportional to 
the slope of the steady-state relation between air flew and engine 
speed and, because these variables increase together, W a ' will 
be positive. 

The first of equations (24) indicates the possibility of negative 
values of t^. Negative values of occur if an increase in en- 
gine fuel flow gives a reduction in the steady-state value of turbine- 
outlet temperature. This possibility exists only near idling engine 
speed when the coefficient bq , which is proportional to the steady- 

state relation between turbine-outlet temperature and engine fuel 
flow, may be negative. In the normal engine operating range, T-^ 

will therefore be positive. 

The relative values of the time constants are of interest, 
because, as previously explained, the ratios of the time constants, 
as compared with unity, determine whether the initial value of the 
indicial response is greater or less than the final value. Infor- 
mation concerning the ratios of the time constants can be obtained 
from a consideration of equations (24). 


The first of equations (24) may be rewritten as 



T 


1-1 V 

1-1 V 


J2 

T 


(27) 


As has been explained, W a * is positive. The coefficient aq, 
which is proportional to the slope of the steady-state relation 
between engine speed and engine fuel flow, is also positive. Equa- 
tion (27) therefore shows that if 


T 


> 


1 



T 


< 1 


then 
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or conversely. If 



Therefore, as shown In the discussion of figure 1 and the 
preceding table. It can be concluded that If a step change In 
engine fuel flow causes the initial value of turbine-outlet tem- 
perature to be greater than the final value, a step change In 
either exhaust-nozzle area or tall -pipe -burner fuel flow will cause 
the init ial value of turbine -outlet temperature to be less than 
the final value. Conversely, if a step change in either exhaust- 
nozzle area or tail-pipe -burner fuel flow causes the initial value 
of turbine -outlet temperature to be greater than the final value, 
a step change in engine fuel flow will cause idle initial value of 
turbine-outlet temperature to be less than the final value. 

The discussion of inai. ai.al response can be Interpreted for 
inputs other than a step input. For example, for the frequency- 
response characteristic, if t^/t > 1 there will be a phase lag 
between input and output and, conversely, if T i/t< 1 there will 
be phase lead between input and output. It therefore follows for 
a turbojet engine with tail -pipe burning for steady-state sinusoidal 
Inputs that, 'if turbine-outlet temperature leads engine fuel flow, 
this temperature will lag either exhaust-nozzle area or tail-pipe - 
burner fuel flow or, conversely, if turbine-outlet temperature 
lags engine fuel flow this temperature will lead either exhaust - 
nozzle area or tail-pipe -burner fuel flow. 

» 

Response far other engine types and other engine variables . - 
As explained in the discussion of equations (3), (4), and (6) in 
the section MALX52S, these equations can be extended to include 
dependent variables other than turbine -outlet temperature and 
engine speed for a turbojet engine with tail-pipe burning, or may 
be extended tp other engine types. 

Equation (27) applies for a turbojet engine with an axial- 
flow compressor and the results regarding indicia! response apply 
only to such an engine. The same results may be obtained from 
equation (6), which holds for an engine with any type of compressor, 
if, from a physical cons iderat ion of the engine, positive or 
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w 

cn 

it* 


negative signs are applied to the coefficients In this equation. 
As has teen explained, the coefficient aq is positive and, in 
general, bj Is also positive, from logic similar to that used 
to determine these signs, it is seen that a^ is positive, b2 • 
is negative, a* is negative, and b_ is positive. With these 
signs applied, equation (6) becomes ° 


*1 

a l 


(V 7 ) 


^2 

®2 




«3 




(28) 


Equation (28) shows that if \/i > 1, then T3/T < 1 or 
t 2A < or conversely. The conclusions previously reached con- 
cerning indioial response and frequency response therefore apply 
to all engines with tail-pipe burning. 

Equation (6) may be applied to engines other than a turbojet 
engine with tail-pipe burning. For example, for the turbine- 
propeller engine, propeller-blade angle can be substituted for 
exhaust -nozzle area in equations (3) and (4). (The tail -pipe -fuel 
flow term vanishes . ) If the coefficients and the time constants 
are redefined, equation (6) can be used to determine whether the 
ratios of the time constants (Tq/r and t 2 /t) are greater or 
less than 1. From a consideration of the steady-state character- 
istics of the turbine -propeller engine, it can be seen that (con- 
sidering corresponding coefficients ) aq, bq, and b2 are posi- 
tive and a2 is negative. With these signs applied, an equation 
similar to equation (28) is obtained for the turbine-propeller 
engine. Therefore, for the turbine-propeller engine, if a step 
change in engine fuel flow causes the initial value of turbine- 
outlet temperature to be greater than the final value, a step 
change in propeller-blade angle will cause the initial value of 
turbine-outlet temperature to be less than the final value, and 
conversely. 

# 

Unpublished engine data show, in general, that for a step 
change in engine fuel flow the initial value of turbine -outlet 
temperature will be greater than the final value and, therefore, 
that for step changes in the other independent variables the ini- 
tial value of turbine-outlet temperature will be less than the 
final value. 


General Control Applications 

The form of the engine transfer functions (equations (3) and 
(4)) can be directly used in setting up any general control 
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configuration. For this purpose, it is necessary that the dynamlo 
characteristics of the engine he known..' Equations (6), (22), and 
(26) can he used to determine muoh of the required information. 
These equations show that -with either of the engine time con- 
stants t or determined from engine dynamic data, the dynamic 
characteristics of idle engine can then he described in -Germs of a 
minimum of steady -state data. These dynamic characteristics can 
then he used in the synthesis of controls. Examples of the pro- 
cedures and the techniques used, in control synthesis are presented 
in references 4 and 5. The utility of- the analysis presented 
herein, however, -will he illustrated hy applying idle analysis to 
scheduled and noninteraction controls. 


Scheduled Controls 


Steady-state conditions for the turbojet engine -with tail- 
pipe burning are determined hy specific settings in the independ- 
ent variables. It follows that, with one variable constant, 
unique relations exist between -the two remaining variables for a 
given steady-state condition. It is on this basis that scheduled 
controls are possible for such an engine. 

Belation between A and Ft* - From, an operational viewpoint, 

it is desirable to operate an engine with tail-pipe burning at maxi- 
mum engine speed and temperature over a range of tail -pipe-burner 
fuel flows and exhaust -nozzle areas . A schedule between exhaust- 
nozzle area and tail-pipe-buraer fuel flow that accomplishes this 
aim is explicitly defined in equations (18) and (23) by the brack- 
eted term, which includes these variables. This relation is 


AA _ ^t_ 

A " 2E 3 F t 


(29) 


As is expected, for steady -state engine operation, exhaust- 
nozzle area and tail-pipe -burner fuel flow act in a similar, but 
opposite, manner. At constant engine fuel flow, the relation 
between exhaust -nozzle area and tail -pipe -burner fuel flow that 
maintains steady-state conditions is given by equation (29) and 
may be used as a basis for a scheduled control. 

Belation between Tg and F 0 . - Ab is expected, a relation 
between turbine-outlet temperature and engine fuel flow exists 
that maintains constant engine speed. This relation is defined in 
equation (19) and is 
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H 2 AT 2 
F e = *e*2 


(30) 


Equation (30) can "be used as a "basis for a scheduled control 
that maintains constant engine speed for a range of exhaust-nozzle 
areas . Such a control is of interest for engine operation when a 
variable-area exhaust nozzle is used as a means of improving the 
thrust -control characteristics of the engine. The schedule in 
equation (30) follows directly from the heat-balance equation for 
the engine and, in a sense, is a restatement of it. 

The schedules presented in equations (29) and (30) oan be 
used in closed-loop control systems in which error correction is 
applied only when the controlled variables deviate small amounts 
from the called -for conditions. 


Noninteraction Controls 

The utility of the analysis will now be illustrated by apply- 
ing the analysis to noninteraction controls. In reference 6, the 
operational form for the engine and control characteristics is 
used both diagrammatically and algebraically to set up a general 
control configuration, to solve for. the conditions required of the 
controllers, to eliminate interaction among the controlled varia- 
bles, and to find the controller operational functions that give 
any desired system response action. 

The problem of satisfactory control of the turbojet engine 
with tail -pipe burning involves the problem of interaction and mul- 
tiple control because there are at least two degrees of freedom, 
or two engine variables, that can be simultaneously controlled. 

This Interaction refers, in general, to the effect of one control 
loop on another control loop.- In reference 6, that type of non- 
interaction is attained whereby any controlled variable setting 
affects only its corresponding controlled variable and no other 
controlled quantity. The method of reference 6 applies to contin- 
uous linear systems and may therefore be applied to the system 
considered herein. 


Development of Control Functions 

From reference 6, a general control configuration is assumed 
and is shown in schematic form in figure 2. The engine and the 


to 
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control arc represented as. matrices , vhere the E's t and C’s are 
engine and control matrix elements , respectively. The input to a 
matrix "box mult iplies every element in Its column. Any output is 
the sum of products in its rev. The control system is generalised 
to this extent: (l) Two dependent engine variables (l and T 2 ) 

are to he controlled and one independent variable (Ft) is bo be 
controlled; (2) the control system employs negative feedback in 
■which errors are applied to the control; and (3) each error is to 
affect every independent variable. 


It is useful to use the following algebraic equations for the 
system of figure 2: 



(31) 


The variables in equations (31) have been placed in a form 
consistent -with the preceding development. The general variables 
in equations (31), in terms of engine variables, are 


AT! 

A X 

A^ 

**e 

T 1 

E 


= *e 

ay 2 

at 2 

ax 2 

AA 

Y 2 

= T 
c 

*2 

* T 



AXg 

AF t 



*3 

"*t 



and 


*11 


R 


21 


a l 

_ ®2 

*13 


TD+1 


" TD+1 

T -jD+1 

- T 2° +1 . 

*23 

T^D+1 

TD+1 1 

3522 ~ 132 

“ TD+1 


r%/x» l 
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It can be seen that the first of equations (31) is simply 
equations (3) and (4) -written in compact form. The second of equa- 
tions f31^ leads to the follotsincr eanations* 



(32) 


The conditions (as. obtained from reference 6) on the controls 
to attain the noninteraction conditions specified are as follows: 


For a setting of H or T 2 to have no effect on Ft, 
tlon (17a) in reference 6 gives 


C 31 " C 32 15 0 


equa- 


for a setting of 
ence 6 gives 


If to affect U only, equation (19) in refer- 
c li E 22 _ ^2 t 2 d+1 

°21 ^1 T i D + 1 


f or a setting of T 2 
ence 6 gives 


to affect Tg only, equation (19) in refer- 

c 12 _ _ Jl2 c H 

c 22 E 11 a l 


and for a setting of Ft to affect Ft only, equation (24) in 
reference 6 and equations (6), (20), and (24) herein give 

C 'l3 S 13 E 22" E 12 E 23 „ 

0 or C'-i* *= 0 

G ' 33 E 11 E 22 ” e 12 E 2 1 

G> 25 _ E 15 E 21" >E H E 25 f3 
C *33 E 11 E 22" E 12 E 21 
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It will be noted, that idle last of equations (33) is the sched- 
ule in equation (29) . For convenience, the control matrix is 
redrawn in figure 3. 

Controlled response aotion . - The preceding noninteraction 
conditions (equations (33),) give the- required ratios between the 
elements of any column of the oonurol matrix (fig. 3) . In order 
to complete the analysis, it remains to choose any one element in 
each column. There is a freedom of independently choosing the 
response of each controlled variable to Its corresponding setting. 
Any controller in the first column will determine engine-speed 
response; any controller in the second column will determine tem- 
perature response ; and any controller in the third column will 
determine tall -pipe-burner fuel-flow response. v The control func- 
tions can be determined from, desired response wherein the response 
functions St are defined as the response of the controlled vari- 
able to its setting. The response functions are therefore 



F-f. 33 ^ 


It follows from reference 6 (equations (36) and (37)) and 
equation (6) presented herein that the control functions in terms 
of the desired response and the engine characteristics are 

c ^22 ^11 ^(TgP+l) 

11 1-^11 E 11 E 22" ]B 12 E 21. 1-%i a l b 2 _a 2' b l 


^22 e 12 ^22 *2 

l-« 2 2 E 12 E 21“ E H E 22 l-# 2 2 a 2 t l" a l‘ b 2 


} (35) 



^*33 

1_ ^ , 33 


J 


CD 

a 



1364 


NACA TBf 2183 


21 


Hhe remaining control functions for noninteraction conditions 
are determined from equations (33) presented herein. It therefore 
is theoretically possible to choose any desired response character- 
istics and to solve for the required control functions to achieve 
this response . 


Analytical results applied to controls . - If the relations 
presented in equations (20) and (24) are combined with the control 
functions as determined from equations (33) and (35), the follow- 
ing equations for "tile control functions result : 


c ll= 7~T“V (TgjEM-X) 
i-^n 

^ii 

p _ ^22 S 2 

C 12“ W~ 

1 “^ f 22 00 

p ^22 H 2‘ a i 

°22= ZT z — 

1-^22 *8 

c , ^*33 S 3 

^33 




1 “^ , 33 


J 


(36) 


Equations (36) give the control functions in terms of the 
desired response and the coefficients, the time constants, the 
steady-state engine data previously discussed. The desired response 
functions ( gt-^i through ^*33) are dictated hy the engine con- 
sol requirements and are chosen by the control designer to fulfi ll 
these requirements. Equations (36) show that if the response func- 
tions are chosen, all of the control functions, with the exception 
°f C21, 'w 111 “be determined from steady-state engine data. The 
, determination of either T-j. or T is required to obtain the con- 
trol function C 21 . 
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Illustrative control . - As an illustration of the method of 
determination of the control function, let the response of engine 

T * 

speed setting n he the response of temperature to 


temperature setting he 


1 and the response of tail -pipe fuel 


l+8gD 


flow to a setting of tail-pipe and fuel flow, he 
for example. 




Therefore, 


^11 m 1 
1 -«n 8 1 D 


The control functions (equations (36)) then "become 


'll 5 


'21* 


1 V J( 1 + ^) 
,_v ££( 1+ ^) 


_ h 2 1 


C 2Z= - 


C’ 23 =- 


H© ag BgD 
*5 1 

ag 8jD 


' 3 3 s 8^ 


7 


(37) 


Equations (37) are the control functions required to give the 
first-order responses chosen. As has been previously explained, 
all of the control functions are determined from steady-state 
engine data, ylth the exception of Ggx* « 


The farm, of the required control functions is given directly 
in equations (37). The control C^x, for example, which is the 
response of engine fuel flow to an engine-speed error (fig. 2), is 
a proportional-plus -integral control with the gain of the propor- 
tional element equal to 
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K 

cn 



and the gain of the integral element equal to 

a 5 t 


The control C 21 , -which is the response of exhaust -nozzle area to 
an error in engine speed, is also a proportional -plus -integral con- 
trol, with the gains on these elements determined as for C^i* It 
is noted that the gain on the integral element of this control can 
"be found from steady-state data, hut that determination of the gain 
on the proportional element requires a knowledge of the engine 
time constant Tq. The remaining control functions (C^ 2 , C 2 g, 
C’ 25 , and C 33 ) are shown to he integral controls with gains as 
determined in equations (37). 


SUMMARY OF RESULTS 

The general form of engine transfer functions for a turbojet 
engine with tail -pipe burning was developed and relations among 
the coefficients and the time constants in these functions were 
found from the transfer functions and from engine thermodynamics . 

1. By use of the developed relations it was shewn that: 

(a) The dominant dynamic characteristics of a turbojet engine 
with tail-pipe burning can he found from steady-state data and one 
transient relation. 

(h) The transfer function that related engine speed to changes 
in engine fuel flow and turbine-outlet temperature is determined 
from steady-state operating data. 

2. The transfer functions, when analyzed to determine indicia! 
response characteristics* showed that, if a step change in engine 
fuel flow causes the initial value of turbine -outlet temperature 
to be greater than the final value, a step change in either 
exhaust-nozzle area or tail-pipe -burner fuel flow will cause the 
initial value of turbine-outlet temperature to be less than the 
final value, and conversely. 
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3. The results of the analysis, -when applied to scheduled 
controls, .gave: 

(a) A relation between exhaust-nozzle area and tail- 
pipe-burner fuel flow that maintains constant engine speed 
and temperature over a range of tail-pipe-bumer operation 

(b) A relation between engine fuel flow and turbine- 
outlet temperature that maintains constant engine speed 
over a range of exhaust-nozzle area 

4. The results when applied to the design of a closed-loop 
noninteraoting control system gave the form of all the control 
functions and a solution in terms of steady-state data 

for six of the seven required control functions. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 29, 1950. 
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APPEJSDIX A 


SIMBGLS 


A 


Ihe following symbols are used, in this report* 
effective exhaust -nozzle area 


a lj a 2* a 3/ . 
' b lt' b 2f' b 3 . 


engine coefficients 


C 


C* 


D 

E 

*e 

Ft 

e 

H 




i 

Kl#K2,K3,K4 

M 


control function to which engine -dependent -variable 
errors are applied 

control function to which engine -independent -variable 
errors are applied 

general coefficient 

differential operator. A- 

dt 

engine-characteristic function 

engine fuel flow in energy units per unit time 

tail-pipe fuel flow in energy units per unit time 

general functions 

acceleration due to gravity 

total enthalpy per pound of air flow 

heat added by engine fuel flow per pound of air flow 

heat added by tail-pipe -burner fuel flow per pound 
of air flow 

polar moment of inertia of engine rotor ' 
coefficients 

Mach number at exhaust-nozzle throat 
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m 

N 

P 

P 

Q 

H 




T 

t 

W a 

X 

T 

a,p 

7 

T j t 1> t 2j T 3 


function of Mach number 
engine rotor speed 
total pressure 
static pressure 
unbalanced engine torque 
gas constant 

response function for controlled engine dependent 
variable 

response function for controlled engine independent 
variable 

total temperature 
static temperature 
a^r flow 

engine independent variable 
engine dependent variable 
general engine time constants 
ratio of specific heats 
control time constants 

engine time constants 


Subscripts : 


1 

2 

3 

e 


compressor inlet (fig. 4) 
turbine exit (fig. 4) 
exhaust -nozzle throat (fig. 4) 
engine 
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f final value 

i initial value 

J,k,v indices 

s 


t 


set value 
tail pipe 



28 


HACA BJ 2183 


APPENDIX B 


are 


DETAILS - OF DEVELOPMENT 
Algebraic Derivation of Beaponse Equations 
Engine speed . - The basic equations used in this development 


Q . jr(TS f p e , A, Ft) 


I M 
dt 


_d_ 

dt 


I ^ (N-N 0 ) 


(1) 

(Bl) 


If Iq is the initial steady -state engine speed, equation (Bl), 
in operational form, is 


Q > I D(AN) 


( 2 ) 


where AN is an incremental change in engine speed. Equation (l), 
expanded and linearized around steady-state operating points, is 


Q AN AF e +^" a AA +^ Ft ' AE t 

where ^" w , for example, is ( ^ ) 

W/fe e ,A, 

bined with equation (2), yields 
ID - ^AN 


(B2) 

. Equation (B2), when corn- 


et 


+ ^.AA + ^.AP. 

e 


e T ~ a““ T " p t *“t 


(B3) 


Equation (B3) can be placed in the form 


3L +1 )" 

r / * 


V + , 

V^7 A 


^ F 


t\ AF i 


*ir »t 


(B4) 
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The solution of the homogeneous equation determines the transient 
response of the engine to any of the forcing functions . The 

term -I /&$ has dimensions of time and Is the time constant of 
the system considered. Thus 


*4 


2 


T «= 



(B5) 


Because equation (B4) Is linear, the principle of superposition 
applies and the response of engine speed to changes in engine fuel 
flow is 


AN 1 f ^*e \ ^e 

N = TD+1 N i F e 

and from equilihrium conditions (D— =*0), it follows that 

*e ^ an/N 

^ * m ^e/ F e 


(B6> 


(B7) 


which is the slope of a steady-state operating curve relating engine 
speed to engine fuel flow at constant A ' and Ft, considered at 
the initial steady-state operating condition. The coefficients of 
the remaining terms of equation (B4) can "be shown to he related to 
the slope of steady-state operating curves in the manner shown hy 
equation (B7). Equation (B4), solved for the response of engine 
speed, then becomes 

AN = a l ^£e 8 2 AA ^ 

B " TD+1 F 0 TD+1 A + TD+1 F t 1 1 

Turbine -outlet temperature or other variables . - The response 
of any other engine variable, such as turbine -outlet temperature, 
will be of a form of the sum of effects due to changes in the inde- 
pendent variables. The response of turbine -outlet temperature to 
changes in the independent variables may be found in the following 
manner; From the hypothesis presented in reference 1, it follows 
that 

(B8) 


Q * ^(N, T2, A, F t ) 
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in a jaanner similar to that used, to oh tain equation (B3) , equa- 
tion (B8) becomes 




y A To 
*2 2 


(B9) 


vhere AA = AE^ = 0. 

Equation (B3), with AA and AEfc equal to zero, becomes 


(■ 


ID - lAN= <F 0 AE e 


(BIO) 


Equations (B9) and (BIO), divided by one another to eliminate 
AB and solved for the response of temperature to changes in engine 
fuel flow, give 


ATo = 




e 


AE, 


e 


(Bll) 


Froip. the development of equation (3), it follows that equation (Bll) 
is of the form 


ATg TjD+1 ^ AF 0 
Tg~ " W" *1 ¥T 


(B12) 


vhere 



and hi is the slope of the steady -state relation between Tg 
and f q at constant A and Ef ®he response of Tg to changes 
in A Anfl F-j; can be found in a similar manner. 

The response of any variable to changes in the independent 
variables can be found in a manner similar to that described for Tg. 
The procedure* for obtaining these expressions can be summarized in 
the following manner: Expressions for unbalanced torque as a func- 
tion of the independent variables are written in which the depend- 
ent variable for which the response is desired is successively sub- 
stituted for each independent variable. These expressions are 
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Q = F e , A, F t ) 

Q = y( TS, T, A, Ft) 
a = *r(H, F e , T, F t ) 
Q = rOS, F e , A, T) 


(B13) 


where T is any dependent variable. 

These expressions for unbalanced torque are then expanded in 
a manner similar to that previously described and placed in a form 
similar to that of equations (B9) and (BIO) . The process used to 
obtain equation (B12) is then repeated. By this procedure the 
response of any variable to changes in the independent variables 
is obtained. 

The responses of engine speed and turbine -outlet temperature 
are 


AN _ a l AF q &2 aa *3 AF t 
* ~ TD+1 F e + TD+1 A + T D+l F t 

ATg T-jD+1 AF 0 TgD+l TgD+l AF-jj 

T 2 TD+1 1 F 0 TD+1 « A TD+1 3 F t 


(3) 

(4) 


where sq through b 3 are slopes of steady-state relations among 
the variables and T through are transient relations. 


Delations amofag coefficients and time constants . - Delations 
among the coefficients and the .time constants in equations ( 3 ) 
and (4) can be obtained by solving these equations for equilibrium 
conditions and for the response of the dependent variables to step 
changes in the independent variables. Equilibrium conditions can 
be found by allowing D to approach zero in equations ( 3 ) nnfl (4). 
These equations then become 


AIlp 

IT 


AT. 


2 1£ 


Tc 


AF, 


ai 


e 


+ a 2 


■e 



+ bg 




AF t 

F t 

Ft 


(B14) 

(B15) 
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ani for step changes In the Independent variables In equations (3) 
and (4), the Initial changes In the dependent variables are obtained 
from equations (3) and (4) by allowing 3) to approach infinity. 



(B16) 


*g 2 ,i 

*2 



(B17) 


An equation similar to equation (B2) can be developed for the 
response of Ig to changes in N, F e , A, and Ft by substi- 
tuting ATg for Q in equation (1). This function expanded and 
linearized around steady-state operating points is 


AT 2 _ K AH j. ^e 


AA 


+ *4 


Aft 

F+ 


(B18) 


Equation (B18) is general and holds far equilibrium conditions and 
for the response of I 2 to step changes in the independent 
variables . 


If equations (B14) and (B15) are combined with equation (B18) 
to eliminate AIT/n and ATg/Tg, the relation for final conditions is 


(bi-aiK].) jT 2, + Cbg-agKi) ~ + (^3" a S K: l) 


AF A a AF t 

K? — 1 + K, M k . _2 
3 A 


(B19) 

aid. in a similar manner, from equations (B16) and (B17 ) , the rela- 
tion for initial response is 



T 


*1 


A£e 

F e 


+ 




*2 




Aft 

*t 


(B20) 

The right sides of equations (B19) and .(B20) are identical and 
therefore the coefficients on similar terms of the left sides are 
equal. If the coefficients on corresponding terms are set equal 
to one another and solved for Kj., the following relations among 
the terms in equations (3) and (4) result; 
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g 




r (T 1“ 1 (T 2’ t) - ~ ( T 3“ t ) (6) 

a l ®2 ®3 


Thermodynamic Equations 


Engine heat 'balance . - Ehe engine-heat-balance equation for 
nonequilibrium operating conditions follows directly by equating 
unbalanced torque times engine speed to the difference "between tur- 
bine and compressor power. If air flow is assumed equal to gas 
flow, the resultant expression is 


E 


e 


*e 

"a 


, TT TT 

~ + h 2 - % 


( 7 ) 


Equation (7) differentiated with H]_ constant is 


dFe Pe dWa 

w a “w a w a 


Q 


dE gE 
"a " w a W a 


+ J- dQ + dff 2 
w a 


(B21) 


For deviation around steady-state conditions, the differentials in 
the equation are differential deviations from steady-state condi- 
tions and the remaining quantities are values of the variables at 
the Initial steady-state condition. If deviations from steady- 
state conditions are considered (Q = 0), equation (B2l) with 
specific heat assumed constant reduces to 


®a 

“» a 


jsr 

p 


dQ 


e 


a 2 ^2 
B e T Z 


( 10 ) 


Equation (10)-, placed in linear form by considering incremental 
changes in the variables, is 


JL A n + 5iL iS. 

»e ” W a "'e B e *2 


(14) 


la this equation, AQ is the difference between final and 
initial unbalanced torque and because initial conditions are steady 
state, AQ is equal to Q. 
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Response of N to F e and T 2 . “ response of N to F e 
and T 2 can toe developed from equation (10) if it is assumed that 
the engine has an axial -flow compressor. The relation between W a 
and for such an engine is 


AW* 

W a 


W a * 


AN 

N 


If this substitution is made in equation (14) and the expres- 
sion for Q given toy equation (2) is also substituted, Q and 
AW a /tf a can toe eliminated from equation (14) to give the following 
expression 


1 + 





%2 AT 2 \ 
He T 2/ 


(B22) 


In equation (B22) the factor IN^/Fq W a ' is the engine time 
constant at constant T 2 and F Q . Because neither A nor Ft 
is involved - in this development, it follows that either of the 
restrictions that Ft is constant or A is constant can toe placed 
in this definition and therefore 


T 


2 “ T 3 


IN 2 

SeWa’ 


( 20 ) 


Equation (B22) solved for the response of N to changes in 
F e or T 2 becomes 


AN 

N 


1 / AF q _ H2 AT^j 1 

T gD+l \F e " H e T 2 / W a ’ 


(19) 


Differentiation and linearization of nozzle equation . - From 
continuity of flow it follows that 


W n 



Logarithmic differentiation of equation (9) with A • 
sidered constant gives ' 


7 3 S 

E 


(9) 


con- 
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dW a dp 3 dA dM 1 dt 3 
+ + 


W e 


P3 A M 2 tg 


(B23) 


■but 


*3 “ M 1 ' 2 


v V 


-1 


which- differentiated is 

d *3 dT 5 _ (r 3 -l) M 2 ^ 
7 3 -l , M 

1 + 


*3 T 3 


also 


(>•¥ 


__3_ 

^3 


Pg * Pg 11 + —S M 

which differentiated and solved for dM^l is 

7 3 -! 


dM jf f*3 _ d P5\ 1 + 2 

M \ P 3 p 3 / 7gM^ 




If the expressions for dtg/tg and dM/M are substituted in 
equation (B23), the following expression results 


dW a = dA , ^3 1 dT 5 t ffgN 1_ / 1_ V 

W a A P 3 " ?T 3 V P 3 ’ ^ y[73VM 2 " ) 


(13) 
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